Low-temperature structure of (^'-Al-Pd-Mn optimized by ab initio methods 
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We have studied and resolved occupancy correlations in the existing average structure model of 
the complex metahic alloy ^'-Al-Pd-Mn [Boudard et al, Phil. Mag. A 74 939 (1996)], which has 
approximately 320 atoms in the unit cell and many fractionally occupied sites. Model variants were 
constructed systematically in a tiling-decoration approach and subjected to simulated annealing by 
use of both density functional theory and molecular dynamics with empirical potentials. To obtain 
a measure for thermodynamic stability, we reproduce the Al-Pd-Mn phase diagram at T=0 K, 
and derive an enthalpy of formation for each structure. Our optimal structure resolves a cloud of 
fractionally occupied sites in pseudo-Mackay clusters. In particular, we demonstrate the presence 
of rotational degrees of freedom of an Alg inner shell, which is caged within two icosahedrally 
symmetric outer shells AI30 and Pdi2. Outside these clusters, the chemical ordering on a chain 
of three nearby sites surprisingly breaks the inversion symmetry of the surrounding structure, and 
couples to an Al/vacancy site nearby. Our refined tiling-decoration model applies to any structure 
within the e-phases family, including the metastable decagonal quasicrystalline phase. 

PACS numbers: 61.44. Br, 61.50.Ah, 61.50.Lt, 71.15.Nc 



I. INTRODUCTION 

The ternary Al-Pd-Mn system is of great interest for 
the study of physical properties of quasicrystals and their 
approximants. Its phase diagram contains both a stable 
icosahedral phase and a stable decagonal phase in the 
vicinity of various other crystalline approximants, many 
of which can be grown as single crystals with high perfec- 
tion. Of particular interest is the series ^^ ^i, • • • of oi- 
thorhombic approximants^ of the icosahedral quasicrys- 
talline alloy z- Al-Pd-Mn, which exists in the Al-rich cor- 
ner of the phase diagram.^ An alternative notation for 
these phases is / = 6, 16, 28, ... , where / denotes the 
strong (0, 0, /) diffraction spot for the interplanar spac- 
ing. ^2 (^23) is frequently named ^-phase. All these 
phases have a common periodicity of 16 A along (0, 1, 0) 
and can be described as ^wo- dimensional tilings perpen- 
dicular to that direction.^ For the ^'-phase (sq) the tiling 
consists of staggered flattened hexagon tiles, where oc- 
casionally phason defects can show up composed of an 
nonagon and pentagon tile.^ There is a similar structure 
where the hexagons are all parallel, denoted ^. In the 
phases, the phason defects are arranged in rows, called 
phason planes and squeezed in between n — 1 rows of 
hexagons forming a layered superstructure. 

The most recent experimental structure refinement of 
was conducted in 1996 by Boudard et al.^ The close 
relationship between ^ and was studied by Klein et 
al. via high resolution electron microscopy.^ So far, no 
structure refinement of the ^^-phases has been reported 
in the literature, but the tiling representation allows to 
extend structure models for ^ and to structure models 
of Important building blocks of all structures are 



pseudo-Mackay icosahedral clusters (PMI) which form 
columns parallel to (0,1,0) and when projected are the 
vertices of the tiling representation discussed above. The 
PMI shares the outer shells (Pdi2 icosahedron and AI30 
icosidodecahedron) with the conventional Mackay icosa- 
hedron (MI) cluster, while the inner shell, an AI12 icosa- 
hedron in the MI, is replaced by a less symmetric shell 
of 8-10 Al atoms. While the positions of the outer 
shells are predicted with high accuracy from diffrac- 
tion data, the positions of the inner Al atoms are less 
well defined and have large uncertainties: in experiment, 
the inner sites are characterized by mixed and partial 
occupancies;^ when projected from a hyper-space model 
they correspond to positions close to the boundary of 
atomic surfaces,^ which suggests that they might be in- 
trinsically unstable. 

Interest in the ^^-phases has spurred when Klein et 
alW discovered a special texture around a partial disloca- 
tion of Burgers vector 1.83 A.^ On the tiling level, ^ the 
dislocation results in the insertion of six phason planes 
of width ^ 171 A, which means it can alternatively be 
viewed as a dislocation in the layered superstructure. 
The dual nature of the dislocation - on the atomic level 
and on the tiling level - and the resulting hierarchy in 
length scales lead to the adoption of the term metadislo- 
cation for this special class of dislocations.'^ Metadisloca- 
tions are known to be important for the plasticity of the 
^^-phases,^ ^ but so far little is known about the atom- 
istic structure of the metadislocation core and the details 
of the metadislocation motion. A refined atomistic struc- 
ture model of the underlying crystal phases as achieved 
in this work can form the basis of a better understanding 
of these two aspects. 
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In this contribution we report a structural refinement 
of the Boudard model of via numerical optimization. 
Simulated annealing with molecular dynamics and em- 
pirical potentials, and density functional theory are used 
to relax the structure into its energy minimum. We de- 
rive an enthalpy of formation for each model and com- 
pare it with the enthalpy of other phases coexisting in 
the Al-Pd-Mn system. This yields a direct measure for 
thermodynamic stability. We identify two main sources 
of disorder that plagued a direct structural determina- 
tion of the structure from the diffraction data. One type 
of disorder is related to the symmetry-broken inner shell 
of the PMI cluster. The symmetry-breaking of this shell, 
along with extra space due to the vacated sites, leads 
to a low-temperature state described by an ensemble of 
configurations, instead of the unique icosahedral one. In- 
dividual representatives of this ensemble are a priori not 
accessible to experimental diffraction refinements, which 
only captures the average structure. We discuss how such 
an average structure may even become ill-defined due to 
missing strong intra-cluster and inter-cluster occupancy 
correlations. The second source of disorder mixes va- 
cancies, Pd and Al atoms in the interstitial parts of the 
structure, not covered by any PMI cluster. 

The symmetry-broken PMI inner-shell configurations 
may acquire special importance due to their potential im- 
pact on the long-range ordering in quasicrystals. For it is 
exactly these atoms that define boundaries of polyhedral 
atomic surfaces defining perfect quasicrystal structures in 
six-dimensional hyperspace.^^ While there is probably no 
hope to accurately determine such individual configura- 
tions in situ in the quasicrystal phase, crystal structures 
like the ^'-phase studied here present a unique opportu- 
nity to study these configurations without the need of 
modelling a quasiperiodic quasicrystalline structure. 

The paper is organized as follows: After a description 
of the initial model and simulation methods in Sec. [Ill 



sites are fully occupied, whereas 15 sites possess occu- 
pancy factors of 0.5 or less. 



and III, respectively, we present details of our results in 
Sec. |IV| first, a low-temperature Al-Pd-Mn phase dia- 
gram is presented in Sec. IV A[ which forms the basis 
for our stability evaluations. The optimization of the 
^'-phase and the rotational degrees of freedom of the in- 
ner PMI shells are discussed in Sec. IIVBI and Sec. IIVCI 
Finally, in Sec.|V|a comparison is made with the isostruc- 
tural AlaPd phase, and first structure models for the pha- 
son defects are presented. We conclude with a short sum- 
mary in Sec. IVl} 



II. 



AVERAGE STRUCTURE OF i'l 
BOUDARD MODEL 



THE 



The space group of the ^'-phase was determined by 
Boudard et al\^ as Pnma (No. 62). The unit cell 
contains about 320 atoms with a composition around 
Al73.5Pd22.4Mn4.i. The lattice parameters are a = 
23.54 A, 6 = 16.57 A and c = 12.34 A. The structure ex- 
hibits a large amount of disorder. Only 19 of 48 atomic 




FIG. 1. (Color online) Atomic clusters in ^^ (a) PMI cluster 
consisting of a centered Mn atom, a symmetry-broken inner 
Al shell, a Pd icosahedron, and an Al icosidodecahedron. (b) 
LBPP cluster composed of two Mn-centered Al icosahedra 
which are capped within a Pd pentagonal prism, and an Al 
icosidodecahedron-type shell. Atomic coordinates and label- 
ing is taken from Ref . 4 Note that in these pictures all atomic 
sites are fully occupied. 

The average structure of consists of two types of 
atomic clusters, the PMI and the so-called large bicapped 
pentagonal prism (LBPP). The clusters are illustrated 
in Fig. [l] The PMI cluster is formed by a symmetry- 
broken (relative to the approximate icosahedral symme- 
try of the outer shell) inner shell, which accommodates 
8-10 Al atoms. Most of the atoms form short (2.35- 
2.5 A), presumably strong bonds with the central Mn 
atom. This inner shell is encaged by the nearly icosahe- 
dral second shell, forming apparently favourable Al-Pd 
patterns with Pd atoms. The second shell is composed 
of a large Pdi2 icosahedron with a radius of 4.44 A, and 
of an AI30 icosidodecahedron with an approximate radius 
of 4.9 A. Both Al and Pd subshells of the second shell 
have only about 0.14 A radial deformation. The signif- 
icance of this cluster is clearly evidenced by the 1.5 A 
gap separating the second shell from the outer atomic 
structure. 

The PMI clusters cover almost 90% of the structure. 
All remaining atoms can be ascribed to a LBPP clus- 
ter. Despite its nearly-perfect pentagonal symmetry, this 
cluster also consists of two nearly spherical shells, that 
are well separated from each other, and also from the 
outer structure. Here, the first-neighbor distances to 
the central atom range from 2.6 to 2.8 A for the inner 
AlioMn2 shell. The outer PdioAl32 shell has about the 
same radii as the second shell of the PMI cluster (4.6 A 
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for Pd atoms, 5.0 A for Al atoms), with a bigger radial 
distortion (about 0.3 A) for the Al subshell. It is still 
well separated from the outer shell, which starts at ap- 
proximately 6.5 A. Only 9 atoms of the LBPP cluster are 
not shared with any adjacent PMI. 




(b) 

FIG. 2. (Color online) Alignment of the atomic clusters in 
^^ (a) Projection of the structure along the stacking axis of 
the PMI clusters. Black lines connect adjacent PMI-cluster 
columns. Dashed rectangle outlines the unit cell, dotted cir- 
cles indicate the positions of the LBPP clusters, (b) Side 
view of one flattened hexagon. For the sake of simplicity, only 
one PMI column is shown. The Pdio pentagonal prisms are 
omitted from the LBPP clusters in this picture. The "free" 
Al atoms not covered by any cluster are encircled by black- 
dashed lines. 

The PMI clusters are stacked on top of each other form- 
ing cluster columns along the b direction. Connecting 
adjacent columns results in a two-dimensional tiling con- 
sisting of flattened hexagons with an edge length of about 
7.6 A. Fig. 0^a) shows a cut through the structure per- 



pendicular to the stacking axis. The LBPP clusters are 
located inside the hexagons and are arranged in a zig-zag 
pattern along the b direction as shown in Fig. [2|b). 

This description comprises more than 98% of the atoms 
in the structure. Only two Al atoms per hexagon (per 
unit cell height) cannot be attributed to any cluster. 

The structure can be divided into three fundamental 
layers perpendicular to the stacking axis: a flat (F) layer 
which is a mirror plane, a puckered (P) layer, and an in- 
verse (I) layer which contains the inversion center. The 
remaining layers in Fig. [2|b) are obtained through sym- 
metry operations of the space group. The subscript i 
denotes inverse images with respect to the I-layer, and 
the asterisk symbolizes mirror images with respect to the 
F- layer. 



III. METHODS 
A. Tiling-Decoration Model 

Our goal was to replace the average structure deter- 
mined from X-ray diffraction with many fractionally oc- 
cupied or mixed sites by an ensemble of configurations, 
each of which is a valid configuration with plausibly low 
energy. Alternatively, recognizing the hierarchical orga- 
nization of the average structure in clusters, that fill the 
space following precisely obeyed linking rules, our ini- 
tial starting models are constructed by a tiling-decoration 
method. The method creates a unique one-to-one corre- 
spondence between tilings and atomic structures. This 
is achieved by the decoration step: each tile, or more 
precisely each predefined tiling object (node, linkage, tile 
interior...) carries an identical atomic motif associated 
with the considered structure. Consequently, we can re- 
fine positional parameters that are associated with each 
decorated tiling objects in a manner analogical to the re- 
finement of Wyckoff positions in a unit cell of a crystal. 
The long-range order of a tiling-decorated structure is 
then defined by the underlying tiling. 

One benefit of this approach is a straightforward use 
of the ^-phase as a replacement of the twice as big 
phase in the refinement process: we simply refine the 
decoration motifs in ^ by total energy calculations, and 
afterwards substitute the underlying tiling representing 
^ by a tiling representing ^^ The tilings for ^ and are 
shown in Fig. [3j In both cases the fundamental tiles are 
Boat tiles (discussed below), packed either in parallel (^) 
or anti-parallel {^') order. 

Decoration rule for the tiling. We associate columns of 
PMI clusters with the "large" nodes of the binary Pen- 
rose tiling (BPT).'^ The zig-zag columns of the LBPP 
clusters project exactly as "small" vertices of the BPT. 
We find that the atomic motifs are in fact consistent with 
a particular subset of rhombic BPT family, namely the 
so-called Hexagon-Boat-Star (HBS) tiling. As shown in 
Fig. [3} the ^- and ^'-phases are pure Boat tilings with no 
Hexagon or Star tiles. The Penrose rhombi, into which 
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I- Layer 



P-Layer 



F-Layer 



FIG. 3. (Color online) Top: HBS tiling models for ^ and 
Red (thick) lines connect adjacent PMI-cluster columns. Unit 
cells are highlighted with dash-dotted lines. Bottom: Tiling- 
decoration rule for the Boat tiles reproducing the Boudard 
model. Larger circles indicate atoms sitting in lower layers. 



the H, B and S tiles can be uniquely decomposed, have 
an edge length of about 6.6 A. We chose a stacking pe- 
riod of 16.5 A for all structures. The initial decoration 
rule for the Boat tiles in the ^- and ^'-tiling was set to 
closely reproduce the atomic structure found by Boudard 
et al We decorate the tiles for each atomic layer sepa- 
rately. It only is necessary to set up decoration rules for 
the three fundamental layers I, P, and F, and align them 
along the stacking axis (I at 0.0, P at 0.125 and F at 
0.25, in units of the unit cell height). The decoration 
rules for the other layers are obtained through symmetry 
operations applied to columnar atomic motifs bound to 
our set of tiling objects. Our decoration rules relating 
the Boudard model of to the Boat tiles are shown in 
Fig. [31 

A generalization of these pure-Boat tiling models to the 
family of HBS tilings requires us to optimize one more 
decoration motif, occuring in the Hexagon tile, when a 
pair of skinny rhombi shares a common edge. In our 
case such a situation leads to a new type of horizontal 
interaction between the LBPP clusters, and must there- 
fore be optimized separately. This final decoration rule 
was optimized (see Sec. VB), and we believe it describes 



faithfully the atomic structures corresponding to the var- 
ious experimentally observed tiling patterns.^ 



B. Total Energy Calculations 

The total energy of a structure is obtained using a 
combination of classical molecular dynamics and density 
functional theory (DFT). First, the candidate structure 
is relaxed into its global energy minimum in a molecu- 
lar dynamics annealing approach. Typically, the sam- 
ple is heated up to 1000 K and slowly cooled down to 
zero Kelvin. We found that an annealing time of about 



50 ps is sufficient for our purposes. The annealing simu- 
lations are performed with the ITAP molecular dynam- 
ics program IMD in the NVT ensemble Atomic inter- 
actions are modeled with the embedded atom method 
(EAM).[I^ Suitable potentials have been developed for 
this work with the force-matching methocff^ using the 
POTFIT program.^ A detailed account of the procedure 
can be found in Ref. |T9, 

The final total energies are calculated with the Vi- 
enna ab initio simulation package VASP . — - - Our cal- 
culations employ the projector- augmented wave (PAW) 
method^^ in the Perdew-Wang generalized gradient ap- 
proximation (PW91-GGA).24 The Brillouin zone IS sam- 
pled with a Monkhorst-Pack grid.^^ The A;-point density 
is chosen in such a way that all energies converge to a 
precision of 10-^ eV. Using a conjugate gradient algo- 
rithm, atomic positions as well as the shape and size of 
the unit cell are relaxed once more, until an accuracy 
of 10~^ eV or better is reached. All calculations are 
performed with a constant plane-wave energy cut-off of 
269 eV. 

To assess a low-temperature Al-Pd-Mn phase dia- 
gram, we calculate total energies of all experimentally 
known stable phases including its binary and pure ele- 
mental subsystems. We define the enthalpy of formation 
AH as the difference of the total energy of a structure 
relative to the composition-weighted total energies of its 
pure elements. Using the program Qhull,^^ a convex 
hull of enthalpy versus composition is calculated. Struc- 
tures that minimize the enthalpy are considered ther- 
modynamically stable (at T=0 K). These phases consti- 
tute the vertices of the convex hull. Edges and facets 
indicate two-phase and three-phase regions, respectively. 
Structures with AH above this "enthalpy surface" are 
unstable. The energy difference AE to the convex hull 
indicates their degree of instability. More details of simi- 
lar VASP and convex hull calculations can be found, for 
instance, in Ref. 28 and 29, 

When comparing unrelated structures, our target accu- 
racy for the total energies is 1 meV/atom, and this accu- 
racy is certainly guaranteed by our convergence criteria. 
In our experience, using for example the Perdew-Burke- 
Ernzerhof generalized gradient approximation (PBE- 
GGA),l^the energy differences AE from PW91-GGA are 
usually reproducible to within 1-3 meV/atom. However, 
when it comes to comparing energy differences for vari- 
ants of the same structure, the energy differences should 
not, in general, be divided by the number of atoms, since 
they only reflect differences in particular places of the 
structure. Such energy differences can be reliably evalu- 
ated with much higher accuracy than 1-3 meV/atom. 
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IV. RESULTS 

A. Low- Temperature Al— Pd— Mn Phase Diagram 

The binary Al-rich Al-Mn and Al-Pd phases are the 
main constraints for the formation of Al-rich ternary 
phases. In the following, we first discuss the binary sub- 
systems as well as all previously known stable ternary 
phases. These phases are summarized in Tab. |T] along 
with the elementary subsystems. Fig. |4] shows the zero- 
temperature "energy-phase diagram" in the Al-rich cor- 
ner without our optimized ^'-phases. 



stable (/IH) 



unstable (JE) 



1 Al (0) 

2 AlnMn (-109) 

3 AlgMn (-181) 

4 AliiMn4(-281) 

5 AlMn (-255) 

6 AbiPdg (-573) 

7 Al3Pd2 (-803) 

8 AlPd (-931) 

9 T-Al-Pd-Mn (-297) 




FIG. 4. (Color online) The Al-Pd-Mn energy-phase diagram 
in the Al-rich corner at T=0 K. Circles label known stable 
binary phases, diamonds indicate known metastable phases, 
squares are either unreported or unknown structures, and tri- 
angles correspond to high-pressure phases. Tie-lines connect 
low-enthalpy phases, constituting the convex hull vertices. 
Enthalpies AH and energy differences AE to the convex hull 
are given in parantheses in meV/atom. The red spot outlines 
the approximate compositions of the ^-phases. 

Al-Pd binary phases. The stable Al-rich phases are 
AUPd, AlsiPdg, AlsPds, and AlPd.^-^^ The latter is a 
FeSi-type structure. The usual 52-type phase (CsCl) 
with the same composition becomes stable only at high 
temperature. Al3Pd2 is an ordered- vacancy state related 
to 52-AlPd, where each Pd atom has only three Pd near- 
est neighbors. The Al2iPd8 phase has a more complex 
structure, and an Al-content similar to the ^'-phase. One 
efficient structural diagnostics for Al-rich Al-TM (tran- 
sition metal) compounds is to look at the local environ- 
ments of the TM atom. These structures often seek to 
maximize the number of Al-TM bonds. TM atoms tend 
to have a few or no TM nearest neighbors. Looking at 
the environments of all TM atoms, (almost) all atoms of 
the structure are represented. Both Pd sites in Al2iPd8 
have a very similar local environment of 10 Al atoms. 



These local environments are reminiscent of the inner Al 
shells in the PMI clusters of the ^'-phase. 

While the A-AUPd phase is reported as a stable Al- 
Pd compound, its structure has never been determined. 
However, the recently determined structure of the AUPt 
was conjectured to be isostructural with it.l^ Substitut- 
ing Pt for Pd in this hexagonal structure, and resolving 
obvious correlations between four fractionally occupied 
Al sites, we find two variants of this phase just slightly 
unstable (by about 5 meV/atom) against decomposition 
to fee Al and Al2iPd8. The structure of AUPd is simi- 
lar to the Al2iPd8: both are packings of similar AlioPd 
clusters. 

At high temperature, the ^^-p hase e xtends into the Al- 
Pd binary system.l ^ti[36i j^^ g^^. V A an optimized low- 
temperature structure model for AlaPd is presented and 
compared with the ^'-phase. 

Al-Mn binary phases. The evaluation of phase sta- 
bilities in the Al-Mn binary phase diagram (based on 
ab initio total energy calculations) results in two sub- 
optimally resolved phases: (i) Al8Mn5 phase: it is diffi- 
cult to refine simultaneously the Al-Mn chemical order- 
ing along with a magnetic structure; (ii) /i-AUMn con- 
tains approximately 570 atoms in a hexagonal unit cell. 
Due to the reported mixed/partial occupancies, an op- 
timization of the low-temperature structure presents a 
significant challenge to the available computational re- 
sources. Thus, our calculated Al-Mn convex hull in- 
cludes four phases: Ali2Mn (AI12W prototype), AleMn, 
AliiMn4, and tetragonal AlMn. ^ Out of these, the low- 
est formation enthalpy AH is achieved by the AliiMn4 
phase, while AlgMn and Ali2Mn enthalpies are barely 
below the tie-line connecting AliiMn4 with fee Al. 

Since the ^'-phase (and all other e-phases) is located 
closer to the Al-Pd binary system in the phase diagram, 
we believe that the Al-Mn binary system has a smaller 
impact on these phases. Even though the AUMn and 
Al8Mn5 phases are refined only suboptimally, the lowest- 
enthalpy structure is AliiMn4. It has also a more similar 
Al content to and should therefore be the most im- 
portant Al-Mn phase to consider. 

Previou sly kn own stable ternary phases include the 
icosahedraP^^ and the decagonal^^ quasicrystalline 
phases, the ^'-phase, and the so-called T-phase, an ap- 
proximant of t he decagonal phase with 12 A stack- 
ing period.^^"^^' The ^'-phase has a composition similar 
to the icosahedral phase, and has been treated as its 
approximant.l^ At the same time it is a proper approxi- 
mant of the decagonal phase with 16 A periodicity. 

The stable icosahedral quasicrystalline phase has a 
composition of Al7i.iPd2o.2Mn8.7.'^ Here we represent 
this phase by the best (lowest energy) approximant model 
available to us, containing 128 atoms in a cubic unit cell 
with an edge length of about 12.6 A, and a composition 
of Al71.9Pd21.9Mn6. 3.^^ The structure model is conven- 
tionally denoted as "1/1". All atoms are incorporated in 
the four PMI clusters and the four Ali2Pd icosahedra, 
decorating the vertices of a canonical cell tiling. The 
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TABLE I. Low-temperature phases in the Al-rich corner of the Al-Pd-Mn phase diagram (cf. Fig. [4]). 
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innermost shells of the PMI clusters in this icosahedral 
approximant are comprised of 10 Al atoms, as opposed 
to our optimal ^'-phases, that contain only 9 Al atoms 
in each inner PMI shell. In our final evaluation of the 
phase stability, this 1/1 approximant is unstable by only 
7 meV/atom. As a competing model, we also computed 
total energies of two models reporte d in l iterature, but 
both had significantly higher energies .^^'^ 

Currently, there is no good estimate of total energies 
for the 12 A decagonal phase. In our database, the struc- 
ture is represented by its approximant T- Al-Pd-Mn. In 
agreement with recent experimental work^^ composition 
and chemical ordering of the optimal low-temperature 
ternary structure is Al3iPd2Mn6.^^ Both space group and 
the number of atoms in the unit cell are identical with 
the reported binary AlaMn high-temperature phase. 

The approximate compositions of the e-phases are 
shown in Fig. |4] Thus, our competing phases are T- Al- 
Pd-Mn, AlsiPdg, and either AlisMn or Al3Pd2. Depend- 
ing on the composition of the candidate structure, we 
will refer the energy difference AE to either tie-triangle. 



Fig. 13 includes the most important low-energy struc- 
tures we optimized in this work. All optimized structures 
are further summarized in Tab. [Ill 

In addition to the reported ternary phases, we also 
include a ternary version of the recently discovered qua- 
ternary structure Al72Pdi8Mn5Si5,^^ with Si substituted 
by Al. According to our stability evaluations at T=0 K, 
the structure is insignificantly more stable than our op- 
timized ^- and ^'-phases, as shown in Fig. [l3) and is 
discussed at the end of Sec. [Vl 



fails. It is well known that quasicrystal-related Al-TM 
compounds often exhibit deep pseudogaps. The posi- 
tion of the Fermi level relative to the pseudogap center 
may have an important impact on the total energy (as 
we indeed show below), and is not dictated necessarily 
by local constraints. At the same time, the occupancy 
factors refined from diffraction data cannot be taken for 
granted, due to the couplings between Debye- Waller and 
occupancy factors. The loose inner shells with low sym- 
metry encapsulated within the highly symmetrical outer 
shells of the PMI clusters may possess nearly continu- 
ous degrees of freedom connecting rather closely spaced 
alternative orientations of this inner cluster. A mix of 
such configurations is nearly impossible to represent by 
a few discrete sites available to a tractable refinement. 



In a second step, we vary the chemical ordering within 
the LBPP clusters. To make the optimization tractable, 
we assumed that the variations of the LBPP cluster espe- 
cially on its 5-fold axis are in the first-order approxima- 
tion independent from the inner-PMI-shell cluster vari- 
ations. As we show later, careful optimization of the 
chemical ordering along this cluster axis revealed strik- 
ing impact on the total energy, and proved that the ^- 
and ^'-phases are stable down to zero temperature, ac- 
cording to our DFT approximated energetics. One of 
the counter-intuitive findings was that the optimal chem- 
ical sequence along the LBPP 5-fold axis places Mn-Pd 
atoms as nearest neighbors - a feature that we have never 
seen in other low-energy structures of the ternary Al-Pd- 
Mn compounds. 



B. Structural Optimization of ^ and ^' with respect 
to T=0 K total energies 

Our structure refinement is divided into two steps. In 
a first step we analyze the occupancies and occupancy 
correlations within the inner PMI shells. Determining 
the correct number of atoms contained within the in- 
ner shells is an important task, in which intuition easily 



Finally, as already mentioned in Sec.|ll| there are some 
atomic sites reported in the diffraction refinement, lo- 
cated in the intersticial spaces between the PMI and 
LBPP clusters. Placement of these atoms is in apparent 
conflict with other Al atoms in the second PMI shells. 
We could never find a plausible low-energy configuration 
that would occupy these sites. 
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1. Occupancies of inner PMI shells 

In this section, we systematically minimize the total 
energy of the ^-phase with respect to (i) the number of 
Al atoms present in the inner PMI shells; (ii) optimal 
shape/orientations of these inner subclusters. The in- 
ner shells are encapsulated within nearly icosahedral sec- 
ond shells, and since they comprise less than 12 atoms, 
their symmetry is necessarily lower than the icosahe- 
dral one, leading to several possible orientations. The 
two constraints selecting among these inner-shell orien- 
tations are vertical (along PMI stacking axis) and hori- 
zontal (perpendicular to PMI stacking axis) interactions 
between the PMI clusters. While the vertical interaction 
is pretty much separable and independent from column 
to column, the horizontal PMI-PMI interactions select 
the favourable spatial patterns of PMI columns. 

The monoclinic unit cell of ^ contains four PMI clus- 
ters (two adjacent PMI columns, each with two clusters 
on top of each other). We start occupying each inner shell 
with 8 Al atoms, and gradually increase the number up 
to 11 Al atoms per PMI. In the following the TM(8) and 
TM(ll) sites of the LBPP clusters are occupied by Mn 
and Al atoms, respectively, consistent with the Boudard 
model [cf. Fig. [ijb)]. The resulting structures contain 
32 Pd and 8 Mn atoms, whereas the number of Al atoms 
varies from 110 for 4x Alg, up to 122 for an inner-shell oc- 
cupation of 4xAlii. The relative stability of each model 
from the appropriate tie-triangle) is evaluated by 
constructing the convex hull from the enthalpies of for- 
mation Ai^ of all binary subsystems as well as all ternary 
phases. 

We found that the structures are most stable with ei- 
ther Alg or Alio shells within each PMI cluster. The 
models are unstable by = 13.9 and 12.6 meV/atom, 
respectively. Slightly higher energies were obtained for 
structures containing both Alg and Alio shells, indicat- 
ing an unfavourable interaction between both shells. The 
PMI clusters with Alg or Alu inner shells are highly un- 
favourable. These structures were unstable by 42.5 and 
33.6 meV/atom, respectively. The Alg inner shells were 
even unstable mechanically, since during the annealing 
process some atoms from the outer AI30 shells drifted in- 
ward, forming Alg shells. On the other hand, the inner 
shells seemed to be "over-occupied" with 11 Al atoms, as 
some atoms drifted toward the outer PMI shells in the 
course of the molecular dynamics annealing. 

Shape of the inner shells. The PMI clusters obtained 
after relaxation are shown in Fig. |5j We observed two 
different types of Alg shells. One type of shell is a nearly 
perfect cube with edge lengths of 2.7-3.0 A. The other 
type of shell can be described as a distorted square an- 
tiprism. The Alg shell is a trigonal prism, capped within 
a larger triangle (tricapped trigonal prism, TTP). Its 3- 
fold axis is aligned (approximately) along the pseudo- 
decagonal (stacking) direction. The Alio shells are best 
described by bicapped square antiprisms (BSA), where 
the 2-fold axes are oriented along the 2-fold axes of the 




4 X Alio 



FIG. 5. (Color online) PMI clusters in ^ after relaxation. The 
vertical axis corresponds to the PMI stacking axis. The outer 
AI30 shells are shown only for the 4xAl8 occupation. The 
triangles of the trigonal AI9 shell and the antiprisms of the 
Alio shells are shaded in blue to outline their 3-fold and 2-fold 
axes, respectively. 



outer PMI shells. The Alg and Alio shells can point in 
different directions as outlined in Fig. |5j For instance, 
in Fig. |5] the 3-fold axis of the TTPs in the upper two 
PMI clusters point along the stacking axis, whereas the 
lower TTPs are aligned along a local 5-fold axis. For 
the Alii shells, no unique description can be established, 
as the shells had variable shape and complicated pattern 
of orientational relationship with respect to the second 
sheh. 



2. Chemical ordering within LBPP clusters 

The inner-shell optimization reported in the previous 
section was performed using a fixed standard model for 
the chemical ordering within the LBPP cluster. In this 
second stage of the refinement, our strategy is to fix the 
energy- minimizing model for the inner shells, and refine 
the LBPP cluster interior. Due to the small energy dif- 
ference between the structures containing Alg and Alio 
shells, we accordingly split the LBPP optimization into 
two branches. Since the Alg variant proved more opti- 
mal, in the following, unless stated otherwise, the PMI 
clusters have Alg inner shells. 

Our optimization first focuses on the atoms on the 
LBPP axis. As already mentioned, the LBPP cluster can 
alternatively be described in terms of separated spheri- 
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(a) (b) (c) 




(d) 

CAI ©Pd •Mn 



FIG. 6. (Color online) (a)-(c) Icosahedra of the LBPP clus- 
ters with different chemical orderings. The structure in (a) is 
unstable by 0.2 meV/atom. The structures in (b) and (c) are 
stable with respect to all competing phases. All three struc- 
tures contain Alg inner shells in their PMI clusters, (d) Align- 
ment of the LBPP cluster from (b) along the PMI stacking 
axis. Only atoms of the flat and adjacent puckered layers are 
displayed in the right picture. Larger circles indicate atoms 
sitting in lower layers. PMI clusters are outlined with dashed 
circles. The Al decagon located on the flat layer is outlined 
with a dash-dotted circle. The unoccupied sites in (b) and (c) 
correspond to Al(21) in Ref. 4 Vac means vacancy. 



cal shells, similar to the PMI cluster. In both clusters 
the second shell consists of two subshells comprised of Al 
and Pd atoms. The inner part is either a pure Alg-n 
shell (PMI), or a composite AlioMn2 shell (LBPP). We 
believe it is more appropriate to change the chemical or- 
dering within the LBPP cluster to obtain a similar atomic 
arrangement as in the PMI cluster. Therefore, we oc- 
cupy the central TM(ll) position of the LBPP cluster 
by Mn and replace the original Mn atoms located at the 
TM(8) sites by Pd as shown in Fig. [6]^ a). In doing so, 
we found a structure lying very close to the convex hull 
{AE = 0.2 meV/atom). The first-neighbor Mn-Pd dis- 
tances within the LBPP clusters are 2.86 A and seem to 
be the key factor for the enhanced stability. The space 
group of this structure is Cmcm (No. 63). The unit cell 
contains 154 atoms in the primitive cell, with a composi- 
tion AlsePdigMua. For comparison, the lowest-energy 
structure we found with the initial chemical ordering 
[see. Fig. [TJb)] was unstable by AE = 12.6 meV/atom. 

We went one step further and broke the mirror symme- 



try of the LBPP cluster by replacing one of the Pd atoms 
within the icosahedra by Al. During the annealing cer- 
tain Al atoms drifted toward the PMI clusters, resulting 
in some Alio inner shells, and leaving one atomic site in 
each LBPP clusters unoccupied as shown in Fig. |6|b). 
The same behavior was observed when both TM(8) sites 
were occupied by Al [see Fig. [6jc)]. The structures were 
almost stable. We removed the diffused atoms in both 
structures to keep the inner PMI shells occupied with 
9 Al atoms, and repeated the annealings. Then both 
structures turned out to be stable with respect to all 
competing phases. Thus, our final optimal models of the 
^-phase contain 152 atoms in the primitive cell, with two 
stable compositions AlsePdiyMna and AlsrPdieMus, re- 
spectively. The space group of both structures is Cmc2i 
(No. 36). 

The electronic density of states of the former struc- 
ture is shown in Fig. j?! Its prominent feature is the 
broad and deep pseudogap. While models with Alg or 
Alio shells have Fermi energies located near the center of 
the pseudogap, as required by the electronic stabilization 
mechanisms, models with unfavorable inner shells, as for 
instance Alg, have Fermi energies shifted by about 0.3 eV 
relative to the pseudogap minimum. 




E-Ep (eV) 



FIG. 7. (Color online) Comparison of the electronic density 
of states (DOS) of ^ with different inner PMI shells. The 
structure with Alg inner shells is stable with respect to the 
competing phases, whereas the structures with Als and Alio 
shells are unstable by 28.0 and 4.7 meV/atom, respectively. 



Finally, we have applied these optimal decoration rules 
to the ^'-tiling, and relaxed the models until forces van- 
ished. The resulting total energies were by 0.4 meV/atom 
lower than their ^-tiling counterparts (for both stable 
compositions). The space group of both structures is 
Pna2i (No. 33). In Ref. [52] we list the relaxed atomic 
positions of our stable ^ and ^' structures. 
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3. Significance of the Al(9) atoms 

So far our models did not include the Al(9) atoms lo- 
cated outside the PMI and LBPP clusters [see Fig. |2] . 
From the diffraction data refinement, the site was half- 
occupied by Al, but its occupancy conflicted with the pair 
of nearby (about 2.3 A) Al(17) atoms, that were reported 
to be fully occupied. Nevertheless, for completeness we 
studied also various models in which this site gets occu- 
pied by an Al atom (while the conflicting Al(17) atoms 
were expected to adjust their position to avoid the short 
distance). 

We observed that during the annealing these atoms 
drifted either into the PMI clusters and increased the oc- 
cupancies of the inner shells, or into the LBPP clusters 
if there were any vacancies available. Only when the an- 
nealing temperature was kept as low as 300-500 K, the Al 
atoms remained at these Al(9) sites. In our stable models 
we fully occupied these sites and relaxed the coordinates 
until forces vanished. The structures were unstable by 
more than 7 meV/atom, corresponding to about 0.5 eV 
cost per Al(9) site. Possibly, these sites are a feature of 
the high-temperature structure. 



C. Internal Degrees of Freedom of Inner PMI 
Shells in ^ and ^' 

The Alg shells forming the core part of the PMI clusters 
have approximate trigonal 3/m symmetry, but they are 
encapsulated within an approximately icosahedral poten- 
tial of the second PMI shells Pdi2Al3o. Due to this bro- 
ken symmetry it is natural to expect that there must 
exist other, energetically nearly equivalent orientations 
of the Alg shell. On the other hand, there are two addi- 
tional strong terms effectively lowering the symmetry of 
the potential inside the PMI: approximately pentagonal 
symmetry breaking due to having one 5-fold axis parallel 
to the global pseudo-decagonal (stacking) axis, and the 
PMI-PMI direct linkages, dictated by the arrangement 
of the PMI-cluster columns. 

The vibrational and configurational degrees of freedom 
related to the reshufflings of the inner-shell atoms are 
activated at finite temperature, and certainly contribute 
to the free energy, providing a stabilization effect. 

Here the dynamics of the inner shells has been studied 
in an ah initio molecular dynamics simulation at 1200 K 
in the NVT ensemble. A time-averaged density plot of 
a 7 A thick slice containing one "layer" of PMI clusters, 
projected along to the PMI stacking axis, is shown in 
Fig. |8] The sharp red-yellow maxima correspond to the 
rigid part of the structure, with atomic positions concen- 
trated around their unique equilibrium positions. The 
weaker, smeared blue-gray areas correspond to extremely 
mobile inner-shell atoms of the PMI clusters. While 
reshuffling, the atoms approximately maintain the trig- 
onal prismatic shape of the shell, resulting in a kind of 
rotational motion: no distinct jumps could be observed. 




FIG. 8. (Color online) Density plot of the atoms in the ^- 
phase at 1200 K, time averaged after 50 ps. The density plot 
is projected along the PMI stacking axis. Dotted lines outline 
the unit cell, dashed circles outline the inner PMI shells. The 
upper part shows a single PMI cluster and its inner Alg shell: 
only atoms within the dashed lines are shown in the density 
plot. For the sake of clarity, one Alg shell is superimposed in 
the density plot. 



To get a first estimate for the energy fluctuations 
as the inner-shell atoms reshuffle, we relaxed the high- 
temperature structure at different times to zero Kelvin. 
The relaxed energies varied up to about 2 meV/atom. 
In all cases the 3-fold axes of the inner Alg shell pointed 
along a local 5-fold axis of the outer PMI shells (not 
necessarily the stacking axis). In addition, we found cer- 
tain configurations where the same inner shell appeared 
slightly rotated with respect to the same 5-fold axis as 
shown in Fig. [9] 



Nevertheless, the time-averaged density plot indicates 
that the preferred orientations of the inner shells must 
be parallel to the stacking axis of the structure. In fact, 
in our lowest energy configurations the 3-fold axes of all 
inner shells point along the stacking direction. 



10 




FIG. 9. (Color online) Snapshots of an inner PMI shell at 
different times after relaxation. The vertical 5-fold axis is the 
stacking axis of the PMI clusters. The triangles of the trigonal 
prisms are highlighted in green to outline their 3-fold axis. 



V. DISCUSSION 
A. ^' in the Binary Al-Pd System 

In 1994 Matsuo et a/P^ reported on an approximant 
to a decagonal quasicrystal with 16 A periodicity, desig- 
nated AlaPd. The orthorhombic unit cell contained ap- 
proximately 280 atoms occupying 300 sites, with lattice 
parameters a = 23.36 A, 6 = 12.32 A and c = 16.59 A, 
and space group Pna2i. The structure is nearly identical 
to the ^^-phase. The main differences are the Mn sites 
that are substituted by Pd, and the inner PMI shells, 
which contain at most 8 Al atoms. A direct comparison 
of the atomic positions can be found in Ref. |4] 

As in the ternary Al-Pd-Mn case, for the purpose of 
optimization of total energies we studied the binary ^- 
phase, rather than the ^'-phase. Interestingly, the inter- 
stitial Al(9) site of the Boudard model was also found 
in AlaPd (Al(21) in Ref. l34|) . The occupancy factor was 
specified as 0.2. Similar to the ternary ^- and ^'-phase, 
we found that these atoms are evidently disfavoured by 
our total energy calculations. 

Varying the number of atoms in the inner PMI shells 
from 8 to 10, we find that like in the ternary phase, the 
optimal shells are Alg with similar trigonal prismatic ar- 
rangements. The second modification that lowered the 
energy by 3 meV/atom was a removal of four Pd atoms 
from four PMI second shells. 

The four PMI clusters contained in the primitive unit 
cell of this binary ^-phase are shown in Fig. [lOj^a). The 
Pd atoms were removed such that the vacancies are lo- 
cated next to each other in the hexagons [cf. Fig.[lO|b)]. 
The centers of the LBPP clusters and the centers of 
the icosahedra within the LBPPs are all occupied by 
Pd: these three Pd atoms form a vertical chain with 
2.86 A distance to each other. The unit cell contains 
112 Al and 38 Pd atoms, corresponding to a composi- 
tion of about Al74.7Pd25.3. At zero Kelvin, the structure 
is by 15 meV/atom unstable to decomposition into fee 
Al and Al2iPd8. At elevated temperatures, instead of 
fee Al, the competing phase should be AUPd. Interest- 
ingly, both of these phases are packings of AlgPd clusters, 
that appear to be a distorted version of the AlgMn clus- 




FIG. 10. (Color online) Optimized AlsPd structure, (a) Side 
view of PMI clusters with central Pd atoms, inner Alg and 
outer Pdii shells, (b) Top view along the stacking axis of PMI 
clusters. Only atoms on the flat and its adjacent puckered 
layers are shown. PMIs are encircled with black-dotted lines. 
Color coding same as in previous figures. 



ters found inside the PMI clusters of the ^-phase. The 
metastability of these binary ^- and ^'-phases is in agree- 
ment with the experimental data. 

In the optimal structure described above, some AlgPd 
clusters are rotated such that their 3-fold axis is no 
longer parallel to the stacking direction of the structure, 
which is a result of molecular dynamics annealing. If 
we insisted they all aligned along the stacking axis, the 
structure is mechanically stable, and its symmetry would 
be consistent with the experimentally determined space 
group Pna2i, but the energy would increase by another 
2 meV/atom. Of course, this agreement should be in- 
terpreted as accidental: even low-symmetry models may 
still belong to a finite-temperature ensemble with higher 
statistical symmetry. 



B. Toward Phason-Line Models 

The e-phase family consists, apart from ^ and of 
a variety of other similar phases. These phases can 
be described by introducing an additional set of tiles, 
a banana-shaped nonagon and a pentagon. As before, 
each vertex is a projection of the PMI clusters along their 
stacking axis. The tiles are always found attached to each 
other, and have been termed phason defect or phason 
line.^ Periodic sequences of these phason lines are called 
phason planes The "larger" e-phases are distinguished 
by the number of hexagon rows that are sandwiched in 
between these phason planes. Their lattice parameters 
along the a and b directions are the same, whereas the 
c parameter depends on the number of inserted hexagon 
rows.^ The simplest of these phases is {sig) which con- 
tains exclusively phason planes with no hexagon tiles. 

Only a few atomistic models exist in the literature 
describing these larger phases. No experimental struc- 
ture refinements have been carried out yet, primarily due 
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to the lack of good quality single crystals. Most inves- 
tigations rely on the average structure of ^^ Tian et 
al. constru cted a ^-phase by means of a cut- and- shift 
method of ^^.'54'55"rpj^^^gj^y^ parts of the same ^' struc- 
ture are translated relative to each other. The model 
is confirmed by comparing the experimental and sim- 
ulated high-resolution transmission electron microscopy 
images. More recently, a structure model for the ^^-phase 
in the Al-Rh system has been proposed using a strong- 
reflections approac/i.^^ Similarly to Tian et a/., the atomic 
positions are deduced from those of the known ^^-phase. 
The space group of both models is B2mm (No. 38). 

Here we use the -phase to obtain a microscopic model 
for the phason lines. We construct the proper HBS tiling, 
and apply our optimized decoration motifs obtained from 
the low-temperature ^- and ^'-phases. In our HBS-tiling 
picture the pentagon is represented by a large Star tile, 
and the nonagon by two Hexagons [see Fig. 11 . Within 





FIG. 11. (Color online) HBS-tiling models for the if, H' and 
-phases. The skinny rhombi are outlined to illustrate the 
different tile-tile interactions existing in each phase. See text 
for more details. 



the tiling the skinny Penrose rhombi are located next to 
each other. This results in a direct contact of the LBPP 
clusters in the structure, a feature that is not present in 
i or i'. 

Before tackling the chain of four skinny rhombi in 
we constructed two other tilings in which the short diago- 
nals of skinny rhombuses form infinite chains [see Fig. 
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Both tilings are composed of Hexagon tiles only. Due to 
their similar arrangement with the hexagons in ^ and 
we termed these tilings H and H' ^ respectively. The 
structure of H consists of PMI planes, formed by connect- 
ing adjacent PMI-cluster columns. The LBPP clusters, 
on the other hand, form further planes which are sand- 
wiched in between the PMI planes. The primitive unit 
cell contains approximately 100 atoms. Similarly, the H'- 
phase is simply an arrangement of staggered PMI planes, 
with the LBPP clusters located in between. The unit cell 
contains twice as many atoms as in the if-phase. Both 
structures have not been observed in the experiment so 
far. Nevertheless, due to their rather small unit cells they 



offer a fast way to calculate total energies for a variety of 
decoration motifs. Afterwards, these motifs can be ap- 
plied to the ^^-phase, and hence to every other e-phase. 

Our lowest-energy structure of the -phase contains 
nearly the same (symmetry-broken) LBPP clusters as our 
optimized ^ and models [cf. Fig.|6|b)]. The inner PMI 
shells are all occupied by 9 Al atoms, forming the usual 
trigonal-prismatic shells. The banana-shaped nonagon 
comprises five entire LBPP clusters per unit cell height. 



as shown in Fig. 12 The nonagon tile, on the other hand, 
contains only one LBPP cluster per unit cell height. The 
structure is about 2.5 meV/atom above our convex hull. 
The primitive unit cell contains 294 Al, 88 Pd, and 16 Mn 
atoms. The space group is Clml (No. 8). The relaxed 
atomic positions of this structure can be found in Ref. [SH 

Our highest-symmetry model for the -phase is unsta- 
ble by 11 meV/atom. The space group is in agreement 
with the experimentally found ones. The structure is 
in principle the same as our best model, except for the 
LBPP clusters: here the Al atoms along the LBPP 5-fold 
axis are substituted by Pd atoms [cf. Fig. |6ja)]. 

For comparison, we also relaxed the structure model 
for the ^-phase, proposed by Tian et al.^ with 540 Al, 
159 Pd, and 36 Mn atoms in the primitive cell. Accord- 
ing to our DFT calculations, their model is unstable by 
more than 30 meV/atom with respect to our competing 
phases. This rather high energy difference is not sur- 
prising: in contrast to our models, their PMI clusters 
contain Alio shells, as well as several (unfavored) inter- 
stitial atoms that are not covered by any cluster. The 
chemical orderings inside their LBPP clusters are also 
in contradiction to our optimized models. Furthermore, 
their banana-shaped nonagons contain only four LBPP 
clusters per unit cell height. 

With the refined models for the flattened hexagons in 
^ and and the optimized models for the phason lines, 
we are now able to construct all variants of the e-phases 
family in the Al-Pd-Mn system. We believe that know- 
ing the optimal decoration motifs for the fundamental 
tiles enables us to construct all other possible structures 
by simply changing the underlying tiling in our tiling- 
decoration model. 

Tab. |TI| summarizes all optimized structures discussed 
in this work by listing compositions, space groups, Pear- 
son symbols and energies. The structures are shown in 
the same order as appearing in Sec. IV B andjVl Our new 
low-temperature energy-phase diagrams with the opti- 



mized structures are shown in Fig. [13] For completeness, 
we also include the recently discovered Al72Pdi8Mn5Si5 
phase in a separate phase diagram, with Si substituted 
by Al yielding a composition of Ali6Pd4Mn. Surprisingly, 
the structure is stable down to T=0 K, with respect to 
all competing phases including our optimized ^'-phases 
[see right part of Fig. 13 . The calculated relative stabil- 



ity of this phase at T=0 K does not change any of the 
conclusions about the optimal structure of the e-phases. 
It should also be remembered that our exploration of the 
low-temperature phase stability remains incomplete due 
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FIG. 12. (Color online) Optimized structure model for ^[ 
with 398 atoms in the primitive unit cell, (a) Projection of 
the structure along the PMI stacking axis. Dashed circled 
outline the locations of the LBPP clusters, (b) Side view of 
one nonagon tile. The Pdio pentagonal prisms and the PMI- 
cluster columns are omitted for the sake of simplicity. The 
nonagon contains five LBPP clusters per unit-cell height. 



tions in chemical ordering and occupancies. The struc- 
ture of these phases is best described as an assembly of 
two types of overlapping clusters, namely the pseudo- 
Mackay icosahedron and the so-called "large bicapped 
pentagonal prism". The fundamental clusters compris- 
ing almost 90% of all atoms are pseudo-Mackay icosa- 
hedra; they differ from the well known Mackay icosahe- 
dron cluster by a reduction of the occupancy and the 
symmetry of the inner shell, from 12 Al atoms (on the 
vertices of an icosahedron) to only 9 Al atoms forming 
trigonal-prismatic shells. For the secondary cluster we 
have determined the chemical ordering along its pentag- 
onal axis which minimizes the energy of the structure. 
We found that the ^- and ^'-phases are, according to our 
calculations, stable down to T=0 K (with a tiny prefer- 
ence for ^'), and might even be a line compound, since 
we found two discrete stable compositions: AlsePdirMua 
and AlsrPdieMna. Both unit cells contain 152 atoms 
(corresponding to 304 in ^'). In an ab initio molecular 
dynamics simulation we showed that the Al atoms of the 
Alg shells inside the PMI clusters are highly mobile which 
can lead to partial occupancy factors as obtained in the 
X-ray refinement. 

A comparison has been made with the nearly isostruc- 
tural AlsPd phase. The structure can be described by 
the same type of clusters. However, in contrast to the 
ternary phases, only metastable structures were found. 

Using a tiling-decoration method we applied our op- 
timized decoration motifs obtained from the refinement 
of the ^- and ^'-phases to a tiling representing the eiQ- 
phase. This phase enabled us to investigate the structure 
of the phason defects in more detail. With these results 
we are now able to construct all variants of the e-Al-Pd- 
Mn phases. 



to the absence of accurate models of the icosahedral Al- 
Pd-Mn phase. 

VI. SUMMARY 

We have systematically optimized the structures of the 
^- and ^'-phases, and established all important correla- 
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FIG. 13. (Color online) Low-temperature energy-phase diagrams of the Al-Pd-Mn system in the Al-rich corner with our 
optimized structures. The structures are labeled according to Tab, [f] and [ill Left: The two ^^-phases are stable down to T—0 K 
with respect to all competing phases. The recently discovered Ali6Pd4Mn phase is not included in the stability evaluation. 
Right: If we include the Ali6Pd4Mn phase in our stability evaluations, both ^^-phases lie slightly above the convex hull. 



TABLE II. Compositions, space groups, Pearson symbols and energies of our optimized structures. The structures are listed 
in the same order as discussed in the paper. A^^oid is the energy difference to our initial convex hull (cf. Fig.[4|, AEne^ is the 
energy difference relative to our new convex hull (cf. left part of Fig. 13). The relaxed atomic positions of structures 16, 17, 
20, 21 and 25 are listed in Ref. 52 
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(eV/atom) 
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^Al55Pdi6Mn4 


Alg 


PI (1) 


aP150 


42.5 


45.8 


-4.7523 


12 


^Al57Pdi6Mn4 


AI9 


PI (2) 


aP154 


13.9 


17.6 


-4.7537 


13 


^-Al59Pdi6Mn4 


Alio 


C12/cl (15) 


mC7316 


12.6 


16.0 


-4.7293 


14 


^Al6iPdi6Mn4 


Alii 


Clcl (9) 


mC324 


33.6 


37.4 


-4.6828 


15 


^AlsePdisMns 


AI9 


Cmcm (63) 


oC308 


0.2 


4.1 


-4.7741 


16 


^AlsePdirMns 


AI9 


Cmc2i (36) 


oC304 


-3.9 


0.4 


-4.7476 


17 


^AlsrPdieMna 


AI9 


Cmc2i (36) 


oC304 


-3.8 


0.4 


-4.7023 


18 


^-Al54Pdi7Mn3 


Alg 


PI (1) 


aP148 


28.0 


31.9 


-4.7426 


19 


^AlssPdirMna 


Alio 


PI (1) 


aP156 


4.7 


8.5 


-4.7133 


20 


^'-AlsePdirMna 


AI9 


Pna2i (33) 


oP304 


-4.3 





-4.7480 


21 


e'-AlsrPdieMna 


AI9 


Pna2i (33) 


oP304 


-4.2 





-4.7027 


22 


^AlsrPdiTMns 


AI9 


Cmc2i (36) 


oC308 


7.0 


10.9 


-4.7238 


23 


^AlssPdieMns 
e'-Al56Pdi9 


AI9 


Cmc2i (36) 


oC308 


7.9 


11.8 


-4.6778 


24 


AI9 


Pna2i (33) 


oP300 


16.6 


16.C1 


-4.5772 


25 


^;-Ali47Pd44Mn8 


AI9 


Clml (8) 


mC796 


2.5 


6.6 


-4.7371 


26 


^;-Ali44Pd47Mn8 


AI9 


Amm2 (38) 


oA796 


11.0 


14.6 


-4.7791 


27 


^-Ali8oPd53Mni2 


Alio 


Amm2 (38) 


0AI47O 


35.0 


38.5 


-4.7429 


28 


Ali6Pd4Mn 


Alio 


Pnma (62) 


0PI68 


-8.5 


-4.7 


-4.6890 



The discovery of a new stable ternary structure does not affect the stability of a binary structure. 
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